A Variational Estimate of the Binding Energy of Charge- Transfer 
Excitons in the Cuprate Superconductors. 
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Abstract 

We present a variational estimate for the binding energy of a Frenkel exciton in the insu- 
lating cuprate superconductors. Starting from the three band Hubbard model we perform 



^ ■ a canonical transformation to 0(t 2 ), where t is the bare nearest neighbour copper-oxygen 



hopping integral. An effective Hamiltonian is then derived to describe the hopping of the 
exciton through the copper oxide plane. The critical parameter in the model is the nearest 
neighbour copper-oxygen coulomb repulsion, V. It is found that a critical value of V is 
needed to observe bound Frenkel excitons, and that these excitons have the same symmetry 
as the parent copper orbital, d x 2_ y 2. We determine the critical value of V using a variational 
approach, and attempt to fit the parameters of the model to known experimental results. 
PACS numbers: 71.35, 74.70V. 



1) Introduction. 



Optical and Raman spectroscopy indicate that charge transfer excitons of the inter- 

atomic kind (d x 2_ y 2 — > p x ,p y ) may be the lowest lying excitations in the insulating cuprate 
superconductors 0]. There are also theoretical predictions that charge transfer fluctuations 
play a crucial role in determining the superconducting and anomalous normal state proper- 
ties of the doped compounds 0-[lO|, although experimental evidence for these excitations in 



the doped phase is less clear. Recent numerical work on the three band Hubbard Model with 
the nearest neighbour repulsion, V, has shown that charge transfer excitations are indeed 



the lowest lying excitations for a critical value of V||ll||. 

We envisage the charge transfer exciton as a localised Frenkel exciton consisting of a 
copper hole excited onto its neighbouring oxygen sites, leaving behind a vacant copper site. 
The resultant bound particle-hole pair then delocalises through the lattice, as illustrated in 
fig. 1. The potential energy of the exciton is the nearest neighbour Coulomb repulsion, V, 
as the oxygen hole has only one neighbouring copper hole. This compares to the potential 
energy 2V for an oxygen hole located amongst occupied copper sites. Thus the condition for 
the exciton to lie in the charge transfer gap is that the potential energy gained in forming 
an exciton must compensate the kinetic energy loss of binding the particle-hole pair. 

In this paper we present a variational estimate for the binding energy of a charge transfer 
exciton in the insulating state of the cuprate superconductors. Starting from the three 
band model, which treats the copper and oxygen orbitals on an equal footing, we perform 
a canonical transformation to 0(t 2 ), where t is the copper-oxygen hybridization integral. 
We thus derive an effective low energy Hamiltonian which describes the charge dynamics 



of the copper-oxide planes. However, since we do not consider terms of 0(t 4 ) we neglect 
superexchange effects. The possible consequences of superexchange are discussed in the 
conclusions. By keeping terms in the Hamiltonian which describe the motion of the exciton 
we derive its energy in both a Neel and a ferromagnetic background. This is therefore a 
variational estimate of the exciton energy. To calculate the exciton binding energy (i.e. 
the energy from the bottom of the conduction band) we compare the exciton energy to the 
minimum energy of the 'single' particle-hole excitation energy. That is, the energy difference 
from the top of the valence band to the bottom of the conduction band. 

The plan of this paper is as follows. In the next section we perform the canonical transfor- 
mation on the three band model. In section 3 we use this to derive the excitonic Hamiltonian 
and calculate the energy of the exciton in the insulating phase. Section 4 discusses the 'single' 
particle Hamiltonians and calculates the free particle-hole gap. In section 5 we calculate the 
condition for obtaining excitons in the charge transfer gap. Finally, in section 6 we attempt 
to fit our model to the experimental results, and conclude. 

2) The Canonical Transformation. 

Our starting point for the expansion is the unperturbed Hamiltonian 



H = H + H u 



(1) 



where 



ija i j 

+ v £ (44^v)- ( 2 ) 

<ij><ja' 



The perturbative part iJ t is defined as 



Ht = ~t £ KW + ^-c), (3) 

<ij>a 

where i and j are copper and oxygen sites respectively, < ij > represents nearest neighbours 
and the operator ^(pj^) creates a copper (oxygen) hole with spin a. A is the charge-transfer 
energy, U d (U p ) is the copper (oxygen) Coulomb repulsion, and V and t are the copper-oxide 
Coulomb repulsion and hybridisation, respectively. 

The effective low energy Hamiltonian is formed using the unitary transformation H = 
e s He~ s , where = —S ||12|| . By expanding in e s it is trivial to show that 

H = H + H t + [S, H ] + [S, H t ] + l[S, [S, H }} + 0(t 3 ). (4) 

To eliminate terms of order t we define S such that H t + [S, H ] = 0. The new Hamiltonian 
is then found by truncating the series at terms greater than 0(t 2 ). Thus H = Hq + H t 2, 
where H t 2 = ^[S, H t ]. H acts on the Hilbert space of \n > where \h >= e s \n > and \n > is 
the basis of our original Hamiltonian. 



By inspection we note that 



i H 0,t (Pjadiv ~ h.C.)} = -t fija(4aPja + h.C.), (5) 

<ij>a <ij>a 

where 

f ija = A+(V-U d )dldi*+(U P -V)pl Pj ,+V ]T p^pjs-V £ <v d .V» ( 6 ) 

<33 i><? 3+3 < n 3>& l 7^« 

and j and j' neighbour i, and i' neighbours j. 
It is relatively straight forward to show that 

S = t J2 (faa)- 1 (p]adia-h.c). (7) 

<ij>cr 

Using the fact that [fija,d ia Pj a ] = [fija,Pj a dia] — H t 2 can be written as 
t 2 

<ij>° <kl>,er' 

- d\*P3<T d { a >Pia' -dlaPjcP^d^ +h.c.\. (8) 

The terms in equation (8) which change the number of oxygen or copper holes are eliminated 
to yield 

t 2 

Ht2= 2 S 51 (/i*0 -1 [pk4r d L'2V - AaPjoPte'dk*' + h 4- ( 9 ) 

<«j>°- <kl>,a' 

Although our Hamiltonian has been simplified it is still very complicated, containing residual 
many body effects resulting from copper-oxygen coulomb repulsion and the effects of onsite 



copper and oxygen coulomb repulsion. 



3) The Excitonic Hamiltonian. 

In constructing the excitonic Hamiltonian we first consider the part of equation (9) which 
hops the oxygen hole around a single copper-oxygen plaquette. This may be written in the 
form 

if = 5>i, (10) 

i 

where 

t 2 

hi = TATV) £ PfoP'*- ( n ) 

v ' <jj'>°- 

Thus when the hole hops to a neighbouring oxygen site via an empty copper site it costs 
an energy r^y\ ■ The symmetry of this Hamiltonian is exploited by noting that the only 
non-zero solution of hi is the constant phase solution P ia = \{p\ c + p\ a + p\ a + p\ a ), where 
the indices on the oxygen operators stand for the four neighbours of the copper site at %. 
From this it can be inferred that our exciton has the same symmetry as the parent copper 
site, d x 2_ y 2, which is the A lg symmetry in Raman notation, hi can be rearranged into the 
more compact form, hi = ,^_ v -. P} c Pi C . The full Hamiltonian is then derived by only keeping 
terms in equation (9) which hop this entity through the 'sea' of Cu 2+ copper sites. We then 
arrive at the final Hamiltonian 



He = E ^ElE ka ,+ ]T (3eie ka + Y,lElE za 

<ik>cra' <ik>a ia 



+ ]T 5ElE lu d\-J k - a + £ eelE ia dld k& + (V + A), (12) 

<ik>cr <ik>a 

where the zero point energy is given as the spin degenerate ground state of one hole per 
copper site. 

The operator E\ a = P} a d i<T creates a spin exciton on site i, while e\ a = P\ a dia spin flips 
and creates a spin 1 exciton on site i. The first two terms in H e hop these excitons through 
the lattice. The remaining terms in H e are self energy terms, which include spin flips. % and 
k are nearest neighbour copper sites. The coefficients are 



t 2 (A - 3V) t 2 V{3V - 5A) + U P (A - W) 

aaa ~ 4A(A + V7 QW ~ 4A{A + V)(A + U p -V) ' 
t 2 _ t 2 (9A-6V) 

' ~ A(A + U p -Vy 7 ~ A{A + V) 

-d S = - e= , A ^ + ^-gy) v ( 13 ) 

4(A + y-c/ d )(A + c/ p -y) 1 ; 



The fact that 5 = — e is due to the effects of exchange. 

With this Hamiltonian we now wish to construct a variational estimate for the energy 
of the exciton. This is done by considering the dispersion of the exciton through either 
a ferromagnetic or a Neel background. In the ferromagnetic case we set a aa > = a aa and 
(3 = 5 = e = 0. The variational wave function is then 



\<t>F >~E e*% n 4tl° >> ( 14 ) 

V-' V cn ; j 



with an energy 



^ <*> = <"' + A) + 3^ry |OTSfe + msh ' ] + Srry ■ (15) 



Similarly, the variational energy of an exciton in the Neel background is found by setting 
a aa > = ol C c and f3 — e — 0. Thus, 



l0iv>=^E e --4,|Neel>, (16) 



cm la 



with 



t 2 (t/ p + y d -2V) t 2 (9A-6VQ 

(A + y-E/dXA + t/p-y) A(A + y) " 1 J 



In the limit of U p — oo E^(k) = E[(k), and for U p < oo E^(min) < E^(min) in our 
allowed parameter range of < V < (A + U p ). It is therefore our intention to consider 
the scenario of an exciton moving in the Neel background, which is also the experimentally 
relevant case. The minimum energy of the exciton can be expressed in terms of the potential 
energy and an effective kinetic energy. For the case of Ud = oo this is written as 



EL = (V + A)-^fZ c (V,U p ,A), (19) 



where 



fzjy, u p ) 



V(9V - 21A) + U P (10A - 9V) + 8A 2 
(A + V)(A + U P -V) 



) 



(20) 



4) Single particle-hole Hamiltonian. 

As we wish to determine the condition for a bound exciton in the charge transfer gap, 
we need to consider the motion of an unbound particle-hole pair in order to find the single 
particle gap energy. In determining the energy of the gap for a free particle hole pair excita- 
tion we intend to study separately the energies of an added hole and a removed hole from the 
insulating phase. The gap energy is then given as E gap = E(N CU + 1)+E(N CU — 1) — 2E(N CU ), 
where N cu is the number of copper sites. We start by considering the energy of a removed 
hole, and then consider the more complicated scenario of an added hole. 

4.1) Removed Hole. 

The Hamiltonian which describes the motion of a removed hole amongst N cu — 1 occupied 
copper sites is]I3] 



H, 



V 



A _ At 2 
2~ ~ A" 



£4 



<ik>a 



Cia + £ 



{A+vy 



At 2 



(21) 



where c\ a creates a particle (destroys a hole i.e. c\ a = d ia ) on site i, and the zero point 
energy is as before. The Nagaoka theorem [|i"4"f informs us that in the strong coupling limit 



the maximum bandwidth for the dispersion of the empty site is obtained in a ferromagnetic 
background. The S = (N cu — 1) branch of the spectrum is then 



A 9t 2 At 2 Rt 2 

E v (k) = - + -^[cosiK) + cos(k y )] - — + _ -. (22) 



The minimum energy of the empty site, at k_ — (n, it), is the top of the valence band 



4.2) Added Hole. 



The Hamiltonian for the added hole follows the analysis of Barford [|15||. In that paper 



the same canonical transformation, equation (4), was performed on the Hamiltonian and the 
effective Hamiltonian was then diagonalised for the i plaquette of one hole on the oxygen 
sites. This was done for the case U p = 0. For a general U p the Hamiltonian is given as 



h c = (2y+|)-(ti-t 2 )E4-^o 

i 

+ ^E^+*iE4-^-Ef^-7x^V (23) 



i,k iS 



A (A + V) 



The zero point energy is again that of a half filled plane (E(N CU ) = — — (A+y-) ) • ^ = it-* 
t 2 = (& + 2v-u d ) anc ^ ^ 3 = A(A+u ) • * represents all the plaquettes with an added hole. The 



operators are defined as 



Si = ^(PlkAl - P k d h) k = 0,+^, vr (24) 



i,fcc 



oE^PL (25) 



2 , 



, i0 = -^=(^4 + 44), (26) 
T i+ = 44 (27) 

and T,_ = (28) 

Sio is the 'Zhang-Rice' singlet operator. The factor ( A+ 2\/-u d ) ai ^ ses fr° m virtual Cu + hop- 
ping, while the remaining terms are accounted for by virtual Cu 3+ hopping. The additional 
complication to the U p = case is the existence of non-zero anti symmetric eigenvalues 
for the singlet Cu 2+ 0~ pair. The oxygen hole can bond in a singlet configuration with its 
neighbouring copper hole and pick up a phase change as it hops around the plaquette. The 
motion of the added hole is complicated because a triplet state on the i site is, in general, 
not orthogonal to a singlet or triplet state on the neighbours of i. It is also unclear how 
the anti symmetric singlet states hybridise with neighbouring sites. It is clear, however, 
that their bands are relatively flat and lie between the singlet and triplet symmetric states. 
Consequently we do not intend to consider these states in our equation. 

We now wish to construct a variational estimate for the minimum energy of this Hamilto- 



nian. To do this we again invoke the Nagaoka theorem|T4| which tells us that the maximum 
bandwidth for the doped hole occurs in a ferromagnetic background. Given this the most 
obvious choice of variational wave function is that of a singlet moving in a ferromagnetic 



background. This has the Bloch wave function 



k> = -7^£ e --4^n4Ti o > 

V^CU I 



J 



This can scatter into a triplet S z — state of the form 



\Pk> = -7^E e --^UTll4Tl >' 



and into a triplet S z — 1 state of the form 



7 fc > = 2-4^= E e--T/ + ^ T 4An4tl°> 

Z V^ CU <m'> j 



These three wave functions give the following overlaps: 



<a k ,\(3 k > = —5 kk >^[cos(k x ) + cos(ky)], 

<a k >\a k > = 6 kk >(l + ^[cos(k x ) +cos(k y )}), 

</3 k .\/3 k > = 5 kk/ {l + ^[cos{k x ) + cos{k y )}) 

and <f3 k i\~i k > = - < a fc /| 7fc0 >= fraclV25 kk '. 



The variational energy is given by E c = <y\ H ^ a k> ^ 

Now, the Hamiltonian acting on the state \a k > gives 



H c \a k >= - (ti - t 2 - t 3 )\a k > -(ti - t 2 - t 3 )^(cosk x + cosk y )\a k 



+ ^-(cosk x + cosk y )\(3 k > — -7=\lk >, (36) 



4^ x 1 —y^- ^ 

so usng the overlaps, equations (32) to (35), E c is given by 



Ec = {2V + A + (-(ti - ta - t 3 ) - {t, - 1 2 - t 3 )f)(l + f ) + ^ + | _ 8 t 2 _ 8t 2 



2 y (1 + f) A (A + V) 

(37) 

where = [cos{k x ) + cos(fc y )]. This value of the energy agrees very well with numerical 



calculations on 4x4 clusters HT3 



In what now follows we will assume that only Cu + virtual excitations are present. This 
is achieved by setting Ud to infinity, i.e. t 2 = 0. Ud is known to be large in the cuprates 
so this assumption has physical validity. The energy of the excitation gap is then simply 
the distance between the top of the valence band and the bottom of the conduction band. 
The top of the valence band is at k = (n, it) and the bottom of the conduction band is 
given at k — (0,0). Thus the minimum gap energy is given as E gap = E c (0,0) + E v (n,7i) 
(E(N CU ) = 0) Although this is not a momentum conserving transistion, we take this as the 
gap energy as we require an upper bound on the value of V for obtaining a bound exciton. 
The gap energy can be rewritten in the form 

E gap = (2V + A) - ^f gap (V, U p , A) (38) 



where 



(39) 



5) Exciton Binding Energy. 



Having calculated the minimum energy of the exciton in the anti ferromagnetic phase, and 
the minimum distance between the top of the valence band and the bottom of the conduction 
band we can now calculate the binding energy of the exciton. This is defined as 

= V + ^(f gap -f exc ). (40) 

The condition E bin = determines the critical value of V for there to be an exciton in the 
charge transfer gap. This is illustrated in fig. 2 where the critical value of V is plotted against 
A (where V = \, A = f and t = ^) for various values of U p . We emphasise that this is an 
■upper bound on the critical value of V, as we have underestimated the gap expected from an 
optical transistion, and the energy of the exciton is variational. Equation (40) illustrates the 
competing effects of the potential energy gain (-V) and the kinetic energy loss ^(f gap — fexc) 
in forming a bound exciton. Notice that the kinetic energy scale is t — ^. In fig. 3 we plot 
the binding energy E bin against V for A = 9t which can be interpreted as a A of 3eV and 
t = l.OeV. 

The effective mass is defined as = A-^f • Thus, the ratio of the effective masses of 
the exciton is given as 

f^ex _ 9 2 E ho i e 8 2 E ex 

m hole ~ 1 dk> ' 1 ' 



i.e. their mass ratio is essentially just the ratio of their band widths. The free hole band 
width is given as E c (tt, tt) — E c (0, 0). The exciton band width is E exc (-n-, n) — E exc (0, 0). This 
gives a mass ratio of 

m ex (AA + 2U P )(A + U P -V)(A + V) 

m hole 2(A + U P )[(V(3V-5A) + U P (A-3V)Y 1 j 



6) Discussions and Conclusions. 

The theory presented in this paper predicts the binding energy of a charge transfer exci- 
ton with d x 2_ y 2 symmetry (Ai 9 in Raman notation) as a function of V, A and U p . We 
now compare these predictions to the experimental data. Recent optical absorption exper- 
iments || give a value of the gap energy as E gap ~ 1.7eV, and a recent Raman scattering 
study of insulating cuprates [Q gives the binding energy of the A\ g exciton as Eu n ~ 0.2eV. 
The value of the effective hopping amplitude t may be estimated from measurements of the 



superexchange interaction (JT6J, PH , J, since J = jp ~ O.lSeV^ and Ud ~ lOeV from pho- 



toemmission data. This gives us an estimate for t ~ 0.6eV. For fixed U p we therefore have 
two equations (19) and (38) with two unknowns and so V and A can be evaluated. This 
proceeds as follows. Firstly t is eliminated from equations (19) and (38) to give the relation 

(E£ c - (V + A))f gap (V, U p ) = (E gap - (2V + A))f e N xc (V, U p ). (43) 



This gives an expression for V and A which can then be reinserted into the expressions for 
the energy of the gap to give the value of t. 

For the case of U p = equation (43) is a quintic polynomial in V and A. This is 
solved using the bisection method to give V ~ 2.7eV and A ~ 3.5eV. This gives the bare 
hopping matrix element as t ~ lAeV. Although these parameters are quite reasonable the 
disappointing feature is that the value of V is large, and greater than U p . This is probably 
unphysical, although direct oxygen-oxygen hopping would be expected to screen U p . It 
should also be remembered that this calculation predicts an upper bound on the value of V 
required for excitons. The mass of the exciton is calculated to be m exc = 0.39m,h o i e . 

For the case U p = oo equation (43) is a quadratic in A which can be solved using the 
usual quadratic formula. This gives a value of V ~ 5.7eV, A ~ OAeV and t = 0.5eV. This 
time all the parameters are inconsistent with the values currently accepted by most to be 
those of the cuprate superconductors, i.e. V ~ leV, A ~ 3eV, t ~ 1.5eV [P|. The exciton 
mass is m exc = 0.36m/j O ; e . 

In conclusion, we have performed a canonical transformation of the three band Hubbard 
model up to 0(t 2 ) in the bare hopping amplitude and derived a variational estimate for the 
energy of an exciton and that of a free particle-hole pair moving through the copper oxide 
plane. We find that the kinetic energies of our particles are rescaled by an effective hopping 
amplitude t — ^r, and that a critical value of the nearest neighbour Coulomb repulsion is 
required to observe bound excitons. This critical value of V is also strongly dependent on 
the value of U p , the oxygen-oxygen Coulomb repulsion. The exciton consists of an oxygen 
hole in a constant phase sum of the orbitals tied to an empty copper site. It therefore 
has d x 2_ y 2 symmetry . Such an excitation is electric dipole forbidden, but an electric- field 



induced dipole transition is allowed if the polarized light is parallel to the electric field JJ. 
Unfortunately such an experiment has not yet been performed. 

In the Raman spectroscopy terminology the exciton has Ai s symmetry and is Raman 
active. These excitations have been observed, and by fitting the experimental data to our 
theory we have found an upper bound for V of 2.7eV assuming U p = 0, or 5.7eV assuming 
Up = oo. The latter result is undoubtedly unrealistic. 

In this paper we have not considered the effects of direct oxygen-oxygen hopping, t pp . 
Including this would have two effects. Firstly it would screen U p . Secondly it would render 
a non-bonding orbital bonding with d 3z 2_ r 2 symmetry. The energy of this exciton is lowered 
while the d x 2_ y 2 symmetry state is raised. An exciton with d 3z 2_ r 2 symmetry is observed in 
both optical and Raman spectroscopy with a binding energy of about 1.7 eV. In addition 
to neglecting t pp we have also neglected the superexchange interaction. This interaction 
would be expected to narrow the valence and conduction bands and hence widen the single 
particle energy gap. However, it is unlikely to effect our estimate of the exciton energy, as we 
assume that it propagates through a Nmboxeel background. Consequently we expect that 
the superexchange term would reduce the value of V required for bound excitons from the 
value calculated by this theory. 
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Figure 1: An exciton in a Nmboxeel background, mving in the copper oxide plane of a 
cuprate superconductor. 

Figure 2: These curves show the boundary between a bound exciton moving in the Neel 
background and a free particle-hole pair for several values of U p . Above the line a bound 
exciton has the lower energy. 

Figure 3: The binding energy of the exciton and gap energy for (a) U p — and (b) 
Up = oo, in units of t, as a function of V. 



